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LONG-TERM  GOALS 

The  long  term  goal  is  regional  modeling  and  data  assimilation,  allowing  us  to  test  models,  and  get 
good  estimates  of  the  state  of  the  ocean.  These  ‘inverse’  estimates  use  the  data  to  overcome  errors  in 
initial  conditions,  in  boundary  conditions,  and  in  dynamical  parameterizations. 

OBJECTIVES 

The  objective  is  to  assimilate  float  or  drifter  data  from  both  coastal  waters  and  the  deep  ocean,  over  a 
period  of  at  least  three  months,  into  a  stratified  primitive  equation  model. 

APPROACH 

Regional  primitive  equation  models  are  ill-posed,  if  the  domain  is  Eulerian.  However,  they  are  well- 
posed  if  the  domain  is  Lagrangian  (  Bennett  and  Chua,1999).  It  is  therefore  possible  to  find  efficient 
solutions  of  the  inverse  or  data  assimilation  problem  in  Lagrangian  domains,  using  the  methods 
developed  here  in  the  last  decade.  My  postdoctoral  fellow  Jodi  Mead  had  committed  the  last  24 
months  to  computing  solutions  of  the  inverse  problem  in  a  Lagrangian  domain  together  with  synthetic 
data.  She  will  continue  work  on  this  project  for  at  least  the  next  12  months. 

WORK  COMPLETED 

Beta-plane  tests  of  the  viscous  shallow  water  equations  in  Lagrangian  coordinates  have  included 
integrations  of  uniformly  rotating  flow  with  open  boundaries,  and  spatially  periodic  integrations  on 
the  time  and  space  scales  of  North  Atlantic  mesoscale  variability. 

The  results  of  the  integrations  are  in  close  agreement  (3%)  with  those  of  Eulerian  integrations  of  an 
otherwise  identical  model,  for  100  days  or  more. 

Lor  variational  assimilation  of  real  N.  Atlantic  float  and  hydrography  data,  the  tangent  linear  form,  and 
adjoint ,  of  the  viscious  Lagrangian  model  has  been  derived  in  spherical  polar  coordinates  .  These 
operators  are  being  implemented  iteratively  in  a  fully  nonlinear  assimilation  scheme. 
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With  help  from  A  Macdonald  of  WHOI,  we  have  used  N.  Atlantic  hydrograhic  data  from 
"HYDROBASE"  to  compute  Montgomery  potentials  on  isopycnal  surfaces.  These  will  provide  'first 
guess'  or  'background'  boundary  conditions  and  initial  conditions  for  assimilations  of  real  float  tracks 
and  depths. 

A  presentation  of  this  work  will  be  given  at  the  meeting  'Lagrangian  Analysis  and  Predictability  of 
Coastal  and  Ocean  Dynamics",  October  2-6,  2000  in  Ischia,  Italy. 

A  rational  derivation  of  open  boundary  conditions  has  been  made  for  diffusive  hydrostatic  and 
nonhydrostatic  dynamics,  as  no  derivation  could  be  found  in  the  published  literature. 

RESULTS 

With  a  uniformly  rotating  flow,  the  viscous  Lagrangian  model  accurately  describes  a  moving  domain 
for  time  periods  of  at  least  160  days.  This  is  a  nested  experiment,  and  the  values  at  the  open 
boundaries  are  interpolated  from  a  known  Eulerian  field  .  Such  interpolation  anticipates  the  use  of 
real  hydrographic  data  provided  on  Eulerian  grids. 

The  integrations  with  time  and  space  scales  of  North  Atlantic  mesoscale  variability  yielded  a  highly 
nonlinear  flow  regime  ,  with  complex  float  trajectories.  Real  float  data  are  of  comparable  complexity. 
Successful  assimilation  of  simulated  data  of  this  type  should  lead  to  successful  assimilation  of  the  real 
data. 

IMPACT/APPLICATIONS 

There  is  little  published  material  on  the  Eagrangian  form  of  the  equations  of  fluid  dynamics,  or 
suitable  numerics,  most  likely  because  the  Eagrangian  viewpoint  is  uncommon.  Semi-Eagrangian 
numerics  are  common  for  efficient  time  integrations  in  numerical  weather  prediction,  but  they  involve 
interpolation  to  a  fixed  grid,  which  is  not  desirable  for  assimilation  of  float  or  drifter  data.  Semi- 
Eagrangian  methods  are  awkward  in  this  context,  but  may  be  computationally  competitive.  The 
Eagrangian  viewpoint  will  give  new  insights  into  ocean  circulation,  both  as  a  forward  model  and  as  a 
basis  for  assimilating  data. 

TRANSITIONS 

Not  applicable 

RELATED  PROJECTS 

We  are  implementing  semi-Lagrangian  methods  to  determine  if  they  conserve  the  well-posedness  of 
the  open  boundaries.  We  are  also  exploring  Lattice  Boltzmann  methods  for  well-posed  regional 
modeling  and  variational  assimilation  in  Eulerian  domains. 
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